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We present a modification to the iodine cell Doppler technique that eliminates 
the need for an observed stellar template spectrum. For a given target star, we 
^ ■ iterate toward a synthetic template spectrum beginning with an existing spec- 
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trum of a similar star. We then perturb the shape of this first-guess template 



cn ' to match the program observation of the target star taken through an iodine 



cell. The elimination of a separate template observation saves valuable telescope 



^ . time, a feature that is ideally suited for the quick-look strategy employed by the 

->^ . "Next 2000 Stars" (N2K) planet search program. Tests using Keck/HIRES spec- 

p . tra indicate that synthetic templates yield a short-term precision of 3 m s~^ and 

a long-term, run-to-run precision of 5 m s~^. We used this new Doppler tech- 
nique to discover three new planets: a 1.50 Mjup planet in a 2.1375 d orbit around 
c§ : HD 86081; a 0.71 Mj^p planet in circular, 26.73 d orbit around HD 224693; and 

a Saturn-mass planet in an 18.179 d orbit around HD 33283. The remarkably 
K^ , short period of HD 86081b bridges the gap between the extremely short-period 

^ ! planets detected in the OGLE survey and the 16 Doppler-detected hot jupiters 
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(P < 15 d), which have an orbital period distribution that piles up at about 
three days. We have acquired photometric observations of two of the planetary 
host stars with the automated photometric telescopes at Fairborn Observatory. 
HD 86081 and HD 224693 both lack detectable brightness variability on their 
radial velocity periods, supporting planetary-reflex motion as the cause of the 
radial velocity variability. HD 86081 shows no evidence of planetary transits in 
spite of a 17.6% transit probability. We have too few photometric observations 
to detect or rule out transits for HD 224693. 

Subject headings: techniques: radial velocities — planetary systems: formation — 
stars: individual (HD 33293, HD 86081, HD 224693) 



1. Introduction 

Of the 171 planets that have been found orbiting stars within 200 pc of the Sun^, the vast 
majority of have been discovered by measuring Doppler reflex motion of their parent stars 
(Butler et al. 2006). Effectively detecting the relatively small Doppler shifts caused by planet- 
sized bodies requires velocity precision of order 1 m s~^. This corresponds to measuring shifts 
in spectral lines to within 10~^ of a resolution element for modern echelle spectrometers such 
as Keck HIRES (Vogt et al. 1994). Such a feat requires intimate knowledge of wavelength 
scale of the spectrometer detector. However, the wavelength scale is not constant with time. 
The instrument is subject to changes on time scales as short as hours or even minutes due to 
effects such as telescope flexure and thermal expansion. It is therefore necessary to measure 
the wavelength scale for each observation using a reference spectrum obtained simultaneously 
with the stellar spectrum. 

One of the most successful and widely-used precision velocity methods is the iodine 
cell technique (Butler et al. 1996, hereafter B96). The B96 technique uses an iodine cell 
placed in front of the entrance slit of an echelle spectrometer, which imprints a dense set of 
narrow molecular lines on each stellar spectrum. Velocity information is extracted by fltting 
a model consisting of a stellar template spectrum multiplied by an iodine reference spectrum 
and convolved with an instrumental proflle, or PSF. The shape of the PSF, wavelength scale 
of the spectrometer, and the Doppler shift of the stellar template are left as free parameters 
in the flt. 

The stellar template spectrum is an observation of the target star obtained without 



^References to published papers and updates on orbital parameters can be found at http://cxoplancts.org 
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the iodine cell, with the instrumental profile deconvolved. Because numerical deconvolu- 
tion algorithms are highly susceptible to noise amplification (e.g. Starck et al. 2002), tem- 
plate observations must have high signal-to-noise (S/R > 500) and high spectral resolution 
(A/ A A > 80,000). Both of these requirements result in longer exposure times compared to 
program observations made through the iodine cell. 

One observing program for which this telescope time overhead is prohibitive is the Next 
2000 Stars (N2K) planet search (Fischer et al. 2005, hereafter Paper I). N2K is an interna- 
tional consortium designed to search for short-period planets (P < 15 d, commonly known as 
"hot jupiters") orbiting metal-rich stars residing within 110 parsecs of the Sun. As telescope 
time is allocated at each of three major telescopes — Keck, Magellan and Subaru — each target 
star is monitored for three consecutive nights. Those stars exhibiting radial velocity (RV) 
variations between 20 and 250 m s~^ are fiagged for follow-up spectroscopic observations 
to search for a full orbital solution. Planet candidates also receive photometric monitoring 
with the automatic photometric telescopes (APT) at Fairborn Observatory (Henry 1999) 
to search for transit events. To date, this "quick-look" observing strategy has resulted in 
the detection of four short-period planets (Paper I; Fischer et al. 2006), one of which, HD 
149026, transits its host star (Sato et al. 2005a). 

Most N2K target stars receive only three RV observations, which means that template 
observations represent more than 25% of the telescope time required for the program. Since 
template observations contain no inherent velocity information, it would be highly advanta- 
geous to bypass them by measuring Doppler shifts with respect to synthetic reference spectra. 
Using this motivation, we have developed a modification to the B96 velocity method whereby 
we derive the shape of the intrinsic stellar spectrum directly from normal iodine cell obser- 
vations. This synthetic, or "morphed," template is then used in place of a deconvolved 
template in the B96 modeling process. 

We present an overview of the forward modeling procedure used to produce morphed 
templates in § 2. In § 3 we show the results of several tests of both the short- and long-term 
precision attainable from morphed templates using Keck/HIRES observations of a set of 
velocity standard stars. In § 4 we detail how this new method is currently being used in 
the N2K planet search. We also report the detection of three new short-period planets, and 
describe the results of our photometric monitoring of the host stars for transit events. Finally, 
in § 5 we summarize our results, propose other potential applications for spectral morphing 
and discuss the new planet discoveries with respect to the growing statistical ensemble of 
short-period exoplanets. 
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2. Spectral Morphing 

The B96 precision velocity technique employs a temperature-controlled pyrex cell con- 
taining gaseous molecular iodine placed at the entrance slit of a spectrometer. The iodine 
lines imprinted on the stellar spectrum serve as a robust wavelength fiducial against which 
the Doppler shift of the stellar lines is measured. The Doppler information from each iodine 
observation is extracted using a forward modeling algorithm described in detail by B96 and 
summarized below. 

Depending on the format of the echelle spectrometer, the stellar spectra are broken into 
600 to 1000 2-A-wide spectral chunks, and each chunk is modeled as 

lobsW = k[Tn{\)DST{\ + AA)] * PSF, (1) 

where /c is a normalization factor, PSF is the model of the spectrometer's instrumental 
response, and "*" denotes a convolution. The fitting procedure is illustrated in Fig. 1. The 
Doppler shift, AA, is applied to the deconvolved stellar template, DST, while the iodine 
reference spectrum, T/2(A), determines the wavelength scale of the spectral chunk. The 
stellar template is multiplied by the iodine template, and the product is convolved with the 
PSF and fit to the observed spectrum. 

The fitting procedure requires two template spectra: the deconvolved stellar template 
and the iodine transmission spectrum. The iodine spectrum has been measured at high 
resolution (R > 200,000) and high signal-to-noise (S/N > 700) using the Fourier Transform 
Spectrometer at Kitt Peak National Observatory. While the same iodine template can be 
used for every RV observation, the stellar template must be measured for each individual 
target star. Stellar template observations are made without the iodine cell in the spectrome- 
ter light path, and are obtained with higher resolution and signal-to-noise ratios than typical 
iodine observations. Templates also require bracketing iodine spectra of rapidly rotating B 
stars, which are used to measure the spectrometer PSF and wavelength scale at the approx- 
imate time of the template observation. The PSF is then removed using a modified Jansson 
deconvolution algorithm (Jansson 1995). 

Our modification to the B96 reduction procedure is to derive the DST directly from an 
iodine observation, thereby circumventing the separate template observation for each target. 
We use a modeling procedure similar to Eqn. 1, replacing the deconvolved stellar template, 
DST, with a synthetic, or "morphed," template spectrum. The shape of the spectrum is 
varied until, when multiplied by the iodine spectrum and convolved with the PSF, a fit 
is obtained to the stellar iodine observation. This morphed deconvolved template is then 
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Fig. 1. — Illustration of the iodine cell Doppler Model. From top to bottom, the shifted 
deconvolved stellar template is multiplied by the iodine reference spectrum. The product is 
then convolved with the spectrometer PSF and fit to the stellar spectral observation obtained 
through the iodine cell. The residuals from the fit are shown at the bottom. 



substituted for the DST in Eqn. 1 during the velocity analysis of the target star's iodine 
observations. 

The fitting procedure requires an accurate initial guess of the shape of the morphed 
DST. A spectrum derived from a model stellar atmosphere would seem like an ideal choice, 
since each individual spectral line could be described analytically. But the fit of Eqn. 1 to 
the observed spectrum is highly sensitive to the exact positions of spectral lines in the DST. 
Small mismatches between model and observed spectra caused by effects such as convective 
blueshift can be compensated by asymmetric adjustments to the PSF model, which can 
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in turn lead to spurious Doppler shifts (Winn et al. 2005). Modern stellar atmosphere 
models cannot accurately account for the identification and exact wavelengths of each of the 
thousands of lines spanning the 1000-A-wide iodine region. 

Instead, we find that an existing observed DST with characteristics similar to the target 
star works the best as the initial guess of the shape of morphed template. Stars with similar 
surface gravity, effective temperature and chemical composition differ only slightly in the 
depth and width of individual spectral lines. We have built a library of DSTs for 42 stars 
spanning the main-sequence and subgiant branch of the HR diagram, with Hipparcos colors 
0.55 < B — V < 1.3 (ESA 1997), and projected equatorial radial velocities Vrot sin i < 
2 km s^^ as measured by Valenti & Fischer (2005). The library templates were observed 
as part of the regular California & Carnegie Planet Search (CCPS, Vogt et al. 2000; Marcy 
et al. 2005b), using the Keck/HIRES echelle spectrometer with the 0.38 arcsecond wide 
E2 decker, which yields a spectral resolution of A/AA ~ 110,000 at 5500 A. Each library 
template observation consists of 3 to 5 sequential exposures, with coadded signal-to-noise 
ratios ranging from 800 to 1000. The PSF is derived from a set of bracketing B star/iodine 
spectra and removed from the template spectrum using the modified Jansson deconvolution 
algorithm (B96). 

The morphing procedure begins by selecting 3 to 6 library template stars with Hipparcos 
colors and absolute visual magnitudes similar to the target star. The non-deconvolved 
template spectra of the candidate library stars are then cross-correlated against the target 
star's "pseudo-template," which is created by dividing the star-plus-iodine observation by 
a B star/iodine spectrum from the same observing run. The library template that yields 
strongest cross-correlation peak, averaged over all orders, is chosen as the initial guess for 
the shape of the deconvolved stellar template. 

With the first-guess spectrum established, the model we use to derive the morphed DST 
is 

/,,,(A) = k[Tj2{X)M{DST)] * PSF (2) 

where M{DST) represents a transformation of the library deconvolved template. We do not 
allow the PSF to vary, but instead derive its shape from a B star/iodine spectrum obtained 
near the time of the observation of the target star. As in Eqn. 1, the DST is shifted by AA. 
But in Eqn. 2 the shape of the spectrum is also allowed to vary by adjusting the depths of 
individual lines and globally broadening the spectrum with a rotational broadening kernel 
corresponding to an equatorial rotation velocity Vrot- 

As with the B96 method, the fit is carried out on individual 2-A-wide spectral chunks. 
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Fig. 2. — Illustration of spectral morpliing algorithm. Top: The first-guess deconvolved 
stellar template (DST, dashed lines) compared to the modified synthetic DST (solid line). 
The synthetic spectrum has been rotationally broadened by 5.1 km s~^and the four spectral 
segments that have been identified as individual lines have been adjusted by a pseudo-optical 
depth factor /j. The borders of the four line segments are marked on the left (dot-dashed 
line) and right (solid). Middle: The DST multiplied by the iodine reference spectrum (solid 
line), convolved by the PSF model, and shifted to fit the stellar observation (dots). Bottom: 
The fit residuals. 



However, in order to modify the line depths of the DST, the spectral chunks are further 
divided into individual line segments using an automated search algorithm. The algorithm 
first locates the minima in the DST that lie further than 4a below the nominal continuum 
level, and then traces along the line wings until either the continuum or a neighboring line 
segment is encountered. Neighboring line segments corresponding to line blends are merged 
and the blend is treated as a single line; a compromise necessary to avoid cross talk between 
the two segments. 

The transformation function can be expressed as 

M{DST) = R{Vrot) * 5^exp{-/,ln[D5T,(A + AA)]} (3) 

where DSTj is the ith spectral segment of the library DST, which is scaled by a factor /j. The 
rotational broadening kernel, R{Vrot), is computed with a linear limb-darkening coefficient 
appropriate for the star's B — V color and the mean wavelength of the iodine spectral region 
(Gray 1992). Fig. 2 shows an example of a morphed deconvolved template, together with 
the resulting fit of Eqn. 2 to an observed stellar iodine observation and the fit residuals. 

Since the library DST is continuum normalized, the logarithmic term in Eqn. 3 can 
be thought of as a pseudo optical depth under optically thin conditions. Our choice of 
this functional form is obviously not motivated by physical considerations since the spectral 
lines in the template are not optically thin. Instead, we chose it because it has the desirable 
behavior of leaving the continuum unchanged while the depth adjustment is applied smoothly 
from the line wings to the core. This has the effect of stretching the line in a somewhat elastic 
fashion. To enforce continuity from one line segment to the next, each segment has a pseudo 
continuum level such that portions of DSTj that lie above this level are left unchanged. For a 
given line segment, the pseudo continuum is the lesser of either the value at the left segment 
boundary or 98% of the stellar continuum. 

In most cases, particularly with N2K targets, we derive the morphed DST from a single 
iodine observation. However, since the fitting procedure is carried out in the rest frame of 
the library star, a refined solution can be obtained by fitting Eqn. 2 to multiple spectral 
observations. The composite morphed DST is formed by replacing the line depths, /», in 
Eqn. 3 with the weighted average. 






T. = ^^4^^, (4) 



where fij is the line depth adjustment for the ith segment of the fit to the jth observation. 
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and aj is the standard deviation of the residuals of the fit of Eqn. 2. The broadening kernel 
is similarly averaged over all observations. 



3. Tests 

The N2K planet search uses a quick-look observing strategy to efficiently search for 
short-period planet candidates within a set of ~ 2000 target stars. Stars are typically 
monitored over three consecutive nights, and those exhibiting conspicuous RV variations — 
typically > 5a — are flagged for follow-up observations. A Jupiter-mass planet in a short- 
period (< 15 d) orbit will produce velocity variations of order hundreds of meters per second, 
necessitating only modest velocity precision. However, the number of extrasolar planets rises 
sharply with decreasing planet mass as dN/dM oc M"^'" (Marcy et al. 2005a), which makes 
it crucial that we attain the highest possible precision in order to detect Saturn-mass planets. 
For example, a 0.3 Mj„p planet in a 15 d orbit around a Solar-mass star induces RV variations 
with an amplitude of only 25 m s~^. Reliably detecting this relatively low-amplitude signal 
from three velocity measurements made during consecutive nights requires a short-term 
measurement precision better than ~ 5 m s~^. 

We are also interested in the long-term performance of morphed templates. While most 
stars are monitored over consecutive nights, the velocity measurements of some stars have 
to be split up over several runs when telescope time is lost to weather . Also, it would be 
ideal if morphed DSTs are stable enough from one run to the next to allow for follow-up 
without having to obtain a traditional template observation. In the following subsections we 
compare both the short- and long-term velocity performance of morphed DSTs to observed 
templates. 



3.1. Standard Stars 

In August 2004, the Keck/HIRES spectrometer was upgraded with a new detector. The 
previous 2K x 2K pixel Tektronix CCD was replaced by an array of three 4k x 2k pixel MIT- 
LL CCDs. The new detector produces much higher velocity precision due to its improved 
charge transfer efficiency and charge diffusion characteristics; smaller resolution elements 
(15 iim pixels compared to the old 24 ixm pixels); higher quantum efficiency; increased 
spectral coverage; and lower read noise. 

Since the detector upgrade, we have monitored a set of 32 chromospherically inactive 
stars as part of the California & Carnegie Planet Search in order to evaluate the RV pre- 
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Fig. 3. — Comparison of velocities derived from a standard observed template (top) and a 
synthetic morphed template (bottom) for HD 185144 (spectral type KOV). The left panels 
show the long-term 64723 velocity precision for each type of template and the right panels 
show the precision during a single observing run denoted by dashed lines in left panels. 



cision of our Doppler reduction pipeline. Each of these RV "standard" stars has received 
intensive monitoring during at least one observing run, with multiple observations spanning 
consecutive nights. Additionally, 14 of the stars have also been observed during at least 
five separate observing runs spanning more than a year and a half. Typical signal-to-noise 
ratios range from 600 to more than 1000, depending on the star's apparent magnitude, with 
a spectral resolution A/AA = 80, 000 at 5500 A. 
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Fig. 4. — The same as Fig. 3 but for HD 9407 (spectral type GOV). 

3.2. Velocity Comparisons 

We used these "standard" stars to perform a differential evaluation of the velocity 
precision of morphed versus observed DSTs. Figure 3 shows an RV plot for one of the more 
intensely observed standard stars, HD 185144 (spectral type KOV). The top panels show the 
velocities measured using an observed DST and the bottom panels show the velocities derived 
from a morphed DST. The velocities measured from a morphed template exhibit a long-term 
RMS scatter of 2.16 m s~^, compared to 1.66 m s^^ for the observed template. The short- 
term RMS for the morphed template velocities is nearly identical to that of the observed 
template (1.44 m s~^ versus 1.56 m s~^). Figure 4 shows a similar velocity comparison for 
HD 9407 (GOV). In this case, the observed-template velocities have RMS < 1 m s~^ for both 
the short- and long-term test cases, while the scatters in the morphed-template velocities 
are about 2 and 3 times larger, respectively. 



Not surprisingly, the mean internal uncertainty, <unc>, of the morphed velocities is 
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Fig. 5. — Histograms comparing the short-term (5 night) velocity precision of morphed 
versus observed (shaded) templates for a set of 32 chromospherically quiet radial velocity 
standards. The mean RMS scatter of the velocities derived from morphed and observed 
templates is given by {o-v,m) ^^^ {^v,m)^ respectively. 



larger than the internal errors produced by the observed templates. This is likely the result 
of small differences between the morphed DST and the true stellar spectrum, which are not 
adequately accounted for by Eqn. 2. These mismatches are also likely responsible for the 
systematic offsets from run to run in the HD 9407 velocity measurements shown in Figure 4 
between days 450 and 550. 

Figure 5 shows histograms of the short-term (single run) velocity scatter, ay, measured 
for the ensemble set of 32 standard stars using morphed and observed templates. Figure 6 
shows a similar comparison for the 14 stars with ~ 1.5 years of observations spanning 
more than five separate observing runs. A comparison of the mean velocity scatter shows 
that morphed templates perform within a factor of 2 of observed templates, with a typical 



-13- 



-i-j 









Long-Term Precision 



<CTv,o> = 2.83 m s" 
<(Tv,M> = 4.23 m s"' 




(Tv [m s ] 



Fig. 6. — The same as Fig. 5 but for a subset of 14 standard stars observed during more 
than five observing runs spanning ~ 1.5 years. 



precision of ~ 3 m s ^ within a single observing run and a long-term precision of ~ 5 m s ^ . 

It should be cautioned that Figures 5 and 6 do not represent unbiased tests of the ab- 
solute precision attainable with the B96 technique. Such a test is complicated by differences 
in the quality of observations from one run to the next, the time coverage and frequency 
of the observations, and the possible presence of undetected low-mass planets. Instead, 
our tests represent a differential comparison of the precision produced by morphed and ob- 
served templates. From this test we can conclude that morphed templates are well-suited to 
the needs of the N2K observing program, providing the necessary short-term precision for 
the initial detection Saturn-mass planets, and long-term precision adequate for monitoring 
planet candidates from one run to the next. 
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4. Application to N2K 

During the 2005B Keck observing semester, we observed 63 metal-ricli stars cliosen 
based on tlie criteria described in Paper I (see also Robinson et al. 2006; Amnions et al. 
2006), using the upgraded HIRES spectrometer. We measured velocities from the spectral 
observations using the B96 Doppler technique, with morphed DSTs created with the first 
iodine observation of each star. 

Figure 7 shows a histogram of the velocity RMS for our sample. For stars with linear 
trends and slopes greater than 10 m s~^ d~^, the scatter is measured with respect to a linear 
fit to the three data points. The majority of the stars (87%) had RMS < 10 m s^'^. The two 
stars with velocity scatter between 10 and 20 m s~^ had predicted velocity "jitter" values in 
excess of 10 m s~^ due to large chromospheric activity indices (Wright 2005). Three stars 
had cry > 250 m s~^, with linear trends indicative of a stellar-mass orbital companion, and 
three stars showed evidence of a short-period planet, with 50 < ay < 250. We obtained 
additional observations of the three planet candidates at Keck Observatory, and present 
here the detection of a 1.5 Mjup planet in a 2.1 d orbit around HD 86081. The remaining 
candidates are still being monitored at Keck. 

Excluding the 8 stars showing excess variability, the mean velocity RMS of our sample 
is 4.47 m s~^. The mean RMS for our N2K targets is larger than the short-term test case 
shown in Fig. 5 due to the lower signal-to- noise of the observations (S/N < 200). 

However, this performance represents a marked improvement over the previous synthetic 
template technique used in Paper I, which used the NSO solar atlas in lieu of our current 
library of DSTs. Our old method produced a mean velocity scatter of ~ 15 m s~^ and 
a long tail in the distribution extending to 60 m s~^. Of the original sample of 211 N2K 
targets monitored at Keck beginning in 2004 and analyzed using our older method, 35 were 
flagged as planet candidates based on having RMS scatter > 3a (Paper 1). These stars 
received additional observations, and we performed the Doppler analysis using observed 
templates to confirm their variability. Three planets from this original sample have been 
reported previously (Paper 1, Fischer et al. 2006). We report here the detection of two new 
planets from our original sample of Keck targets: a 0.72 Mjup planet orbiting HD 224693 
with circular 26.7 d period and a Saturn-mass planet in an eccentric 18.2 d orbit around 
HD 33283. 

In the following sections we present the stellar parameters of the host stars, orbital 
parameters of the planets, and the result of our search for transit events. 
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Fig. 7. — Distribution of RMS velocities measured for 64 N2K target stars. Three binaries 
(SBl) and three planet candidates with 20 < ay < 250 m s~^ lie outside of the plot range. 
The three stars with ay > 9 m. s^^ have jitter estimates in excess of 10 m s^^. With 
these 8 outliers removed, the mean is (o"y) = 4.47 m s"^, with 87 % of the stars having 
(T\/ < 10 m s"^. 

4.1. HD 86081 

HD 86081 is an F8V star with V = 8.73, B — V = 0.664, a Hipparcos parallax-based 
distance of 91 parsecs (ESA 1997), and an absolute visual magnitude. My = 3.93. The star 
is chromospherically inactive, with no emission seen in the core of the Ca II H & K lines. 
We measure Shk = 0.16 and logR'^j^ = —5.03 and derive a rotational period of 40.1 d 
following the calibration of Noyes et al. (1984). Our high resolution spectroscopic analysis, 
described in Valenti & Fischer (2005), yields T^ff = 6028 K, log^ = 4.36, [Fe/H] = +0.257 
and Vrot siiii = 4.2 km s~^. Interpolation using the Yi et al. (2001, , hereafter Y^) isochrones 
provided a stellar mass estimate of 1.21 ±0.05 Mq and a stellar radius of 1.22 ±0.1 Rq. The 
stellar parameters are summarized in Table 1. 



We started monitoring HD 86081 in 2005 November, and our Doppler analysis using 
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Fig. 8. — Top: Phased radial velocity measurements of HD 86081 (observed template). 
The dashed line shows the best-fit orbital solution. Bottom: The residuals of the orbital fit 
(filled circles) and the difference between the velocities measured with observed and morphed 
templates (open triangles). The RMS of the fit residuals is 3.20 m s^^ and the RMS difference 
between the two sets of RV measurements is 2.8 m -sT^ . 



a morphed template yielded a non-linear velocity variations with RMS = 154 m s~^ from 
the first three measurements. Additional monitoring revealed a clear, 2.1375 d periodicity. 
We list our velocity measurements in Table 2 along with the Julian Dates and internal 
measurement uncertainties. We report here only the velocities measured using a high-quality, 
observed template, due to the higher velocity precision. 

In order to search for a the best-fit orbital solution, we augmented the internal mea- 
surements errors by adding in quadrature a photometric jitter estimate of 3.7 m s^^ based 
on the chromospheric activity index of the star (Wright 2005). Figure 8 shows the phased 
RV observations of HD 86081 together with the best-fit orbit solution: P = 2.1375 d, 
K = 207.7 km s~^, and e = 0.008 ± 0.004. From the orbital parameters and adopted stellar 
mass of 1.21 Mq, we derive a minimum planet mass Mpsini = 1.50 Mjup and semimajor 
axis a = 0.039 AU. The lower panel of Fig. 8 shows the fit residuals (filled circles), which 
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Fig. 9. — Phased radial velocity mesurements of HD 224693 and best-fit orbital solution. 
The fit residuals have an RMS of 4.07 m s~^. 



have an RMS scatter of 3.20 m s"^ and yxl = 0.81. Also shown in the lower panel are the 
differences between the measurements made with morphed and observed templates (open 
triangles), which show excellent agreement (RMS = 2.8 m s~^). 

We list the best-fit orbital parameters in Table 5 along with the parameters determined 
from the morphed-template velocities. The two sets of parameters are identical within errors. 
The parameter uncertainties were estimated using a Monte Carlo method. For each of 1000 
trials, the fit was subtracted from the measurements, the scrambled residuals were added 
back to the measurements, and a new solution was determined. The standard deviations of 
the parameters derived from all trials were adopted as the la uncertainties. 



4.2. HD 224693 

HD 224693 is aV = 8.23, G2IV star with B -V = 0.639, a parallax-based distance 
of 94 pc and My = 3.36 (ESA 1997). The star is chromospherically quiet, with logi?'^^ = 
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—5.15 and a rotational period of 27.4 d, based on the calibration by Noyes et al. (1984). Our 
spectral synthesis modeling determines Tg// = 6037 K, \ogg = 4.38, Vj-ot sini = 3.5 km s^^and 
[Fe/H]= +0.343. The mass of the star derived from Y^ isochrones Yi et al. (2001) is estimated 
to be 1.33 ± 0.1 Mq with a radius of 1.70 ± 0.3 Rq. 

HD 224693 was observed at Keck as part of the N2K program beginning in 2004 July. 
The first four Doppler observations showed a modest velocity RMS of 16 m s~^. At that time, 
the synthetic template algorithm was still being tested, so a standard template was made 
to confirm the low amplitude variation. We obtained a total of 24 radial velocities, listed in 
Table 3. We augmented the internal uncertainties with a jitter estimate of 4.0 m s~^ based 
on its chromospheric activity index (Wright 2005). 

The velocities are well fit by a Keplerian with P = 26.73 d, K = 40.2 m s~^, and 
eccentricity, e = 0.05 ± 0.03. Figure 9 shows the phased radial velocity measurements along 
with the best-fit Keplerian model. The RMS fit to the model is 4.07 m s~^ with y^xE^ = 0.85, 
consistent with the measurement errors. From our orbital parameters and the adopted stellar 
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mass of 1.33 Mq we derive a minimum planet mass Mp sin? = 0.71 Mjup and semimajor 
axis a = 0.233 AU. The full list of derived orbital parameters is given in Table 6. 



4.3. HD 33283 

HD 33283 is a \/ = 8.05, G3V star with B -V = 0.641. Based on the Hipparcos 
parallax, it has a distance of 86 pc and absolute visual magnitude My = 3.36 (ESA 1997). 
The star is chromospherically inactive with S = 0.10, logi?^^ = —5.60 and a predicted 
rotational period of 55.5 d. Spectral synthesis modeling yields Te// = 5995 K, logs' = 4.21, 
V^otsini = 3.2 km s~^, [Fe/H] = +0.366. Interpolation of the Y^ (Yi et al. 2001) isochrones 
suggests a mass of 1.24 ± 0.1 Mq and radius of 1.20 ± 0.1 Rq. We summarize the stellar 
characteristics in Table 1. 

We began monitoring HD 33283 at Keck in 2004 July. The first five radial velocities 
showed a modest RMS of 21 m s~^, with a 17 m s^^ change between the first and second 
observation. At that time, our synthetic template algorithm was still being tested, so a 
standard observed template was obtained to confirm the variations. The 25 velocities mea- 
sured with the observed template are listed in Table 4 along with the internal measurement 
uncertainties. Based on the chromospheric emission index, we expect 4.0 m s~^ of stellar 
jitter, which we added in quadrature to the internal uncertainties. 

The best-fit Keplerian orbit solution yields K = 25.2 m s^^, P = 18.179 d, and e = 
0.48. The phased radial velocities and orbit solution are show in Figure 10. The fit has 
RMS = 3.60 m s^^ and ^/xi = 0.77, consistent with the measurement errors. From our 
adopted stellar mass and best-fit orbit parameters, we estimate a minimum planet mass of 
0.33 Mjup and a = 0.168 AU. The full orbit solution is summarized in Table 6. 



4.4. Photometry of the Host Stars 

Between 2005 November and 2006 February, we obtained high-precision photometry of 
two of the three planetary host stars in this paper with the T8 and TIO 0.8 m automatic 
photometric telescopes (APTs) at Fairborn Observatory. The APTs can detect short-term, 
low-amplitude brightness variability in the stars due to rotational modulation in the visibility 
of magnetic surface features such as spots and plages (e.g., Henry et al. 1995) as well as 
longer-term variations associated with stellar magnetic cycles (Henry 1999). The photometric 
observations can help to establish whether observed radial velocity variations are caused by 
stellar activity or planetary-refiex motion (e.g., Henry et al. 2000a). Queloz et al. (2001) 
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and Paulson et al. (2004) have found several examples of periodic radial velocity variations 
in solar-type stars caused by photospheric spots and plages. The APT observations are also 
useful to search for possible transits of the planetary companions (e.g., Henry et al. 2000b; 
Sato et al. 2005a). 

The T8 and TIO APTs are both equipped with two-channel precision photometers each 
employing two EMI 9124QB bi-alkali photomultiplier tubes to make simultaneous measure- 
ments in the Stromgren b and y passbands. The APTs measure the difference in brightness 
between a program star and a nearby comparison star. The local comparison stars used for 
the program stars were HD 84664 and HD 203 for HD 86081 and HD 224693, respectively. 
Stromgren b and y differential magnitudes were computed and corrected for differential ex- 
tinction with nightly extinction coefficients and transformed to the Stromgren system with 
yearly mean transformation coefficients. The external precision of the differential magni- 
tudes is typically between 0.0012 and 0.0017 mag for these telescopes, as determined from 
observations of pairs of constant stars. Finally, we combined the Stromgren b and y differen- 
tial magnitudes into a single {b + y)/2 passband to maximize the precision of the photometric 
measurements. The individual photometric observations of HD 86081 and HD 224693 are 
given in Tables 7A and 7B, respectively. Further information on the automatic telescopes, 
photometers, observing procedures, and data reduction techniques can be found in Henry 
(1999) and Eaton et al. (2003). 

The photometric results for the two stars we measured are summarized in Table 8. 
We have no photometry as yet on HD 33283. Columns 2 and 3 give the date range of 
the photometric observations and the total number of individual differential magnitudes, 
respectively. The standard deviations in column 4 refer to the spread of the (6 -|- y)/2 
measurements around the mean values of the data sets. All standard deviations are consistent 
with the constancy of the host stars and their comparison stars. Periodogram analyses of the 
two data sets do not reveal any significant periodicities. We computed the semi-amplitudes 
of the light curves (column 5) with least-squares sine fits of the observations phased to the 
radial velocity periods. The resulting amplitudes are both 0.0004 mag or below; this absence 
of light variability on the radial velocity periods supports planetary-refiex motion as the 
cause of the radial velocity variations in both stars. Although we have no photometry of 
HD 33283, it is the most inactive of the three stars (Table 1) and so should show a similar 
lack of photometric variability. 

The sixth column of Table 8 gives the geometric probability of transits for the three 
systems we observed photometrically, computed from equation 1 of Seagroves et al. (2003) 
and assuming random orbital inclinations. The final column gives the results of our search 
for planetary transits of the two companions as shown in Figures 11 and 12 and discussed 
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Fig. 11. — Top: The complete set of Stromgren {b + y)/2 photometric observations of 
HD 86081 obtained with the T8 0.8 m APT at Fairborn Observatory and plotted against 
orbital phase of the planetary companion. The predicted time, depth, and duration of 
possible transits are shown schematically. The star exhibits no optical variability on the 
radial velocity period larger than 0.0004 mag or so. Bottom: The observations around the 
predicted time of transit are replotted with an expanded scale on the abscissa. The error 
bars are described in the text. Transits deeper than 0.001 mag or so are ruled out by these 
observations. 
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Fig. 12. — Top: The complete set of Stromgren {b + y)/2 photometric observations of 
HD 224693 obtained with the TIO 0.8 m APT at Fairborn Observatory and plotted against 
orbital phase of the planetary companion. The predicted time, depth, and duration of 
possible transits are shown schematically. The star exhibits no optical variability on the 
radial velocity period larger than 0.0003 mag or so. Bottom: The observations around the 
predicted time of transit are replotted with an expanded scale on the abscissa. The error 
bars are described in the text. Our observations are too few to rule out possible transits. 
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below. 

The top panel in each of Figures 11 and 12 plots the complete photometric data set for 
each star against the planetary orbital phase, computed from the orbital elements in Tables 5 
and 6; zero phase in each case refers to a time of inferior conjunction (mid-transit). The 
solid curve in each panel approximates the predicted transit light curve, assuming a planetary 
orbital inclination of 90° (central transits). The out-of-transit light level corresponds to the 
mean brightness of the observations. The transit durations are calculated from the orbital 
elements, while the transit depths are derived from the estimated stellar radii (Table 1) and 
planetary radii computed with the models of Bodenheimer et al. (2003). The bottom panel 
of each figure shows the observations around the predicted time of mid transit replotted 
with an expanded scale on the abscissa. The horizontal error bar below each predicted 
transit curve represents the approximate ±lcr uncertainty in the time of mid transit, based 
on Monte Carlo simulations and the uncertainties in the orbital elements. The vertical error 
bar represents the precision of the measurements in each case. 

Our photometric observations of HD 86081 cover the predicted transit window quite 
well. The 21 observations within the transit window in the bottom panel of Figure 11 have 
a mean brightness level of 1.3003 ± 0.0004 mag. The 114 observations outside the transit 
window have a mean of 1.2993 ± 0.0002 mag. Thus, the two values agree to within 0.001 
mag, or ±2cr, placing very shallow limits on possible transits. In the case of HD 224693, our 
data around the predicted time of transit are far to sparse to rule out possible transits at 
this time (Figure 12) and we are continuing to monitor this star. 



5. Summary and Discussion 

We have presented a modification to the B96 precision radial velocity technique that 
eliminates the need to obtain template observations for each target star. We replace observed 
templates with synthetic versions derived by modifying the deconvolved stellar template 
(DST) of a star with similar characteristics as the target star. The depths and widths of 
the spectral lines in the surrogate DST are modified until, when multiplied by an iodine 
template and convolved with a model of the spectrometer PSF, a fit to the target star's 
observed star-plus-iodine spectrum is obtained. Based on our tests of Keck/HIRES spectra 
of a set of chromospherically inactive stars, we have determined that morphed templates yield 
velocity measurements with precision within a factor of two of those produced by observed 
templates; with a short-term precision of ~ 3 m s~^ and long-term stability within 5 m s~^. 

While observed templates yield the highest velocity precision, the observational time 
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savings afforded by morphed DSTs make them favorable alternatives for many types of planet 
search programs. Examples include surveys containing evolved stars such as K giants, where 
stellar jitter is the dominant noise source (Frink et al. 2002; Sato et al. 2005b); surveys 
conducted on small telescopes where it is not possible to obtain high enough signal-to-noise 
for template observations; testing for binarity among SIM grid star candidates (Frink et al. 
2001; Gould 2001); and surveys containing faint stars such as M dwarfs. 

Morphed templates are also ideal for the "Next 2000 Stars" (N2K) planet search, which 
uses a quick-look observing strategy to efficiently search for short-period planets within a 
large sample of stars (Paper I). Synthetic templates are an integral part of the N2K observing 
strategy and have already yielded important results. We report here the discovery of three 
new planets from the Keck N2K program: a 1.50 Mjup planet in a 2.1375 d orbit around 
HD 86081; a Saturn-mass planet in an eccentric, 18.179 d orbit around HD 33283; and a 
Jovian-mass planet with a circular, 26.73 d orbit around HD 224693. 

The quick-look strategy used by the N2K consortium involves observing stars during 
three consecutive nights, and Monte Carlo simulations have shown that the strategy is most 
sensitive to planets with Mp sini > 0.5 Mjup and periods between 1.2 and 14 days (Paper I). 
However, the planets orbiting HD 33283 and HD 224693 have orbital periods significantly 
longer than 14 days. These detections were aided by the fact that intensive follow-up 
observations were not initiated until after the fourth or fifth observations were obtained. 
These two cases serve as important illustrations of how the sensitivity of the N2K quick-look 
observing strategy can be greatly enhanced with additional observations and an increased 
time baseline. 

The 2.1375 d orbital period of HD 86081b bridges the gap between the 3 day pile- 
up seen in the period distribution of hot jupiters (HJ, 3 < P < 15 d; Butler et al. 2006; 
Udry et al. 2003) discovered by RV surveys, and the extremely short period (P < 2 d) 
planets found by the Optical Gravitational Lensing Experiment planet transit search (OGLE; 
Konacki et al. 2003; Bouchy et al. 2004). While these "very hot jupiters" (VHJ) are common 
amongst the OGLE detections, comprising 3 out of 5 transit detections, Doppler surveys have 
previously found only two other Jovian-mass planets with a period significantly less than 3 
days (HD 73256, P = 2.55 d and HD 189733, P = 2.22 d Udry et al. 2003; Bouchy et al. 
2005). Gaudi et al. (2005, hereafter G05) showed that this apparent inconsistency between 
RV and transit searches can be reconciled by considering the different selection biases for the 
two search methods and the P^^/^ sensitivity scaling for ground-based transit searches. By 
accounting for these effects, G05 concluded that VHJs represent a rare class of exoplanet, 
with an occurrence rate only 10% that of HJs, based on the aggregate sample from RV and 
transit searches. 
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Since RV surveys such as the Cahfornia & Carnegie Planet Search (CCPS) and N2K are 
complete down to Saturn-mass planets with orbital periods less than 15 days {K = 25 m s~^), 
we can use the combined sample of planets from these two surveys to estimate the relative 
number of VHJs and HJs. Out of ~ 1800 total stars, these two RV surveys have yielded 
16 HJs and two VHJs (using a somewhat arbitrary lower cutoff at P = 3 d and therefore 
including HD 83443 as done by G05). Thus, we estimate the ratio of HJs to VHJs is ~ 8, 
which is in quantitative agreement with the estimates of G05. A much more detailed analysis 
will be provided in a future publication in the N2K Consortium series (Fischer & Laughlin 
2006). 

The relative paucity of VHJs provides an interesting puzzle, and the answer is likely 
related to the way in which close-in giant planets are formed and how they interacted with 
their early circumstellar environments. One of the leading theories suggests that giant planets 
migrate inward along nearly circular orbits through interactions with their circumstellar 
disks. These interactions can take the form of gravitational scattering and physical collisions 
with planetesimals in the remnant disk (Murray et al. 1998), or through type I or II orbital 
migration (Trilling et al. 1998; Ida & Lin 2004). These disk interaction scenarios predict very 
few high-mass planets {Mp > 2 Mjup) at small semi-major axes due to mass-loss suffered 
from Roche Lobe overflow. This prediction is in accordance with the observed deficit of 
high-mass, short-period planets. 

There is another, very different theory for forming short-period giant planets. In this al- 
ternate scenario, planets with initially high eccentricity and small periastron passages become 
tidally circularized at small orbital separations. Using the observed properties of extrasolar 
planets. Ford & Rasio (2006) showed that there exists a sharp cutoff in the mass-period 
distribution nearly equal to twice the Roche limit, corresponding to the ideal circularization 
distance (Faber et al. 2005; Gu et al. 2003; Rasio & Ford 1996). If the original eccentricity is 
induced by planet-planet or planet-star interactions (e.g. Weidenschilling & Marzari 1996; 
Rasio & Ford 1996), then the expectation is that there would be fewer high-mass hot jupiters 
due to their inertial resistance to scattering, similar to the migration scenario and also in 
agreement with the observed data. 

The Roche limit scales inversely with planet-star mass ratio and proportionally with 
the planet radius. Since the planet radius is nearly constant for Mp > 1 Mjup, the tidal 
circularization theory helps explain the rarity of massive planets in extremely short-period 
orbits. Indeed, the two most recently discovered VHJs, HD 189733b and HD 86081b, lie 
safely beyond the ideal circularization radius. However, Ford & Rasio (2006) caution that 
the distance of the cutoff was derived under the assumption that the slope follows the Roche 
limit since the slope cannot be constrained by the current observational data. Additional 
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leverage can be gained by increasing the number of targets monitored in order to find rare 
high-mass planets close to the circularization limit, or increasing precision in order to find 
low-mass, short-period systems. This highlights the importance of programs like N2K that 
are specifically designed to efficiently search for short-period planets within a large sample 
of stars. 
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Table 1. Stellar Parameters 



Parameter 


HD 86081 


HD 224693 


HD 33283 


V 


8.73 


8.23 


8.05 


Mv 


3.93 


3.36 


3.36 


B-V 


0.664 


0.639 


0.641 


Spectral Type 


F8V 


G2IV 


G3V 


Distance (pc) 


91 


94 


86 


Fe/H 


0.257 (0.03) 


0.343 (0.03) 


0.366 (0.03) 


Teff (K) 


6028 (44) 


6037 (44) 


5995 (44) 


Vrot sin i km s^^ 


4.2 (0.5) 


3.5 (0.5) 


3.2 (0.5) 


\ogg 


4.36 (0.07) 


4.38 (0.07) 


4.21 (0.07) 


M, (Mo) 


1.21 (0.05) 


1.33 (0.1) 


1.24 (0.1) 


-R* {Rq) 


1.22 (0.1) 


1.70 (0.3) 


1.20 (0.1) 


Shk 


0.16 


0.13 


0.10 


logi?W 


-5.03 


-5.15 


-5.60 


Prot (d) 


40.1 


27.4 


55.5 
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Table 2. Radial Velocities for HD 86081 



JD 


RV 


Uncertainty 


-2440000 


(m s-i) 


(m s-i) 


13694.156 


-168.30 


4.20 


13695.155 


129.92 


4.29 


13697.155 


52.46 


4.77 


13746.991 


151.71 


3.40 


13747.997 


-184.74 


3.20 


13748.936 


201.56 


2.72 


13749.891 


-209.26 


3.06 


13750.907 


189.65 


2.12 


13750.936 


195.23 


1.85 


13751.883 


-170.66 


2.09 


13751.937 


-196.19 


3.68 


13752.948 


165.28 


2.30 


13752.986 


172.90 


2.22 


13753.060 


197.83 


2.11 


13753.124 


201.11 


2.25 


13753.912 


-121.22 


2.83 


13753.934 


-132.36 


2.27 


13753.962 


-142.86 


1.97 


13754.033 


-171.86 


2.17 


13754.069 


-185.19 


2.07 


13775.897 


-142.57 


3.60 


13776.835 


180.54 


3.50 


13777.920 


-191.64 


3.46 


13778.927 


193.64 


3.43 


13779.925 


-207.33 


3.80 


13781.064 


190.15 


3.64 
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Table 3. Radial Velocities for HD 224693 



JD 


RV 


Uncertainty 


-2440000 


(m s-i) 


(m s-i) 


13191.097 


47.16 


5.09 


13198.131 


20.71 


5.54 


13199.091 


17.74 


5.50 


13200.042 


9.35 


5.06 


13367.705 


-37.79 


2.65 


13368.713 


-26.67 


3.26 


13369.714 


-26.54 


2.50 


13370.736 


-13.68 


11.3 


13550.112 


-20.60 


3.48 


13551.122 


-18.08 


3.25 


13552.090 


-27.37 


3.38 


13571.088 


32.34 


3.31 


13603.089 


-15.32 


3.06 


13604.062 


-21.39 


2.87 


13605.021 


-24.05 


2.87 


13692.906 


-6.92 


1.97 


13693.850 


-3.39 


1.95 


13694.762 


3.36 


2.11 


13695.757 


9.33 


2.09 


13696.751 


23.32 


2.01 


13723.810 


31.15 


2.65 


13724.756 


35.05 


2.68 


13746.696 


-6.67 


3.14 


13752.761 


40.52 


2.67 
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Table 4. Radial Velocities for HD 33283 



JD 


RV 


Uncertainty 


-2440000 


(m s-i) 


(m s-i) 


13014.852 


10.57 


5.45 


13015.858 


1.88 


5.35 


13016.850 


-11.57 


4.98 


13071.767 


-29.19 


4.85 


13072.843 


-31.26 


5.07 


13368.955 


4.10 


4.94 


13369.808 


4.87 


4.99 


13397.796 


-8.86 


4.52 


13398.870 


-20.06 


5.14 


13428.792 


17.14 


4.50 


13605.139 


10.66 


4.35 


13692.917 


-12.34 


3.93 


13693.979 


-2.42 


4.23 


13694.975 


-0.03 


4.35 


13695.994 


6.55 


4.33 


13696.981 


9.04 


4.26 


13723.925 


9.44 


4.49 


13724.932 


4.19 


4.53 


13746.852 


-15.10 


4.81 


13747.891 


-2.94 


4.60 


13748.809 


-0.13 


4.37 


13749.824 


5.66 


4.27 


13750.822 


12.91 


4.14 


13751.875 


14.61 


3.34 


13752.878 


13.17 


3.89 
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Table 5. Orbital Parameters for HD 86081 From Observed and Morphed Templates 



Parameter 


Observed 


Morphed 


P(d) 


2.1375 (0.0002) 


2.1378 (0.0002) 


V (JD) 


2453694.8 (0.3) 


2453694.4 (0.4) 


e 


0.008 (0.004) 


0.007 (0.005) 


Ki (m s-i) 


207.7 (0.8) 


207.8 (1.0) 


uj (deg) 


251 (40) 


192 (70) 


Mp sin i {Mjup) 


1.50 


1.50 


asini (AU) 


0.039 


0.039 


Fit RMS (m s-^) 


3.20 


4.38 


V^ 


0.81 


0.87 


Nob, 


26 


26 



^Time of periastron passage. 
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Table 6. Orbital Parameters for HD 224693 and HD 33283 



Parameter 


HD 224693 


HD 33283 


P(d) 


26.73 (0.02) 


18.179 (0.007) 


T/ (JD) 


2453193.9 (3.) 


2453017.6 (0.3) 


e 


0.05 (0.03) 


0.48 (0.05) 


Ki (m s-i) 


40.2 (2.) 


25.2 (2.) 


uj (deg) 


6. (200) 


155.8 (8.) 


Mp sin i {Mjup) 


0.71 


0.33 


asini (AU) 


0.233 


0.168 


Fit RMS (m s^^) 


4.07 


3.60 


V^ 


0.85 


0.77 


Nobs 


24 


25 



^Time of periastron passage. 
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Table 7. Photometric Observations of HD 86081 



Observation Date 


A(6 + y)/2 


(HJD - 2,400,000) 


(mag) 


53,739.9097 


1.2972 


53,739.9168 


1.2987 


53,739.9654 


1.2985 


53,740.8213 


1.2976 


53,740.9088 


1.2964 



Note. — Table 7A, along with Ta- 
ble 7B (for HD 224693), is presented 
in its entirety in the electronic edi- 
tion of the Astrophysical Journal. A 
portion is shown here for guidance 
regarding its form and content. 
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Table 8. Photometric Results for the Planetary Host Stars 





Date Range 




a 


Star 


(HJD - 2,450,000) 


Nobs 


(mag 


(1) 


(2) 


(3) 


(4) 



Semi-Amplitude Transit Probability 

(mag) (%) Transits 

(5) (6) (7) 



HD 86081 3739-3777 135 0.0019 0.0004 ± 0.0002 17.4 no 

HD 224693 3691-3733 22 0.0018 0.0003 ± 0.0004 4.2 ? 



